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We study the distribution of the Schmidt coefficients of the reduced density matrix of a quantum 
system in a pure state. By applying general methods of statistical mechanics, we introduce a 
fictitious temperature and a partition function and translate the problem in terms of the distribution 
of the eigenvalues of random matrices. We investigate the appearance of two phase transitions, one 
at a positive temperature, associated to very entangled states, and one at a negative temperature, 
signalling the appearance of a significant factorization in the many-body wave function. We also 
focus on the presence of metastable states (related to 2-D quantum gravity) and study the finite 
size corrections to the saddle point solution. 
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I. INTRODUCTION 

Entanglement is an important resource in quantum in- 
formation processing and quantum enabled technologies 
[U |2]. Besides its important applications in relatively 
simple systems, that can be described in terms of a few ef- 
fective quantum variables, it is also widely investigated in 
many-body systems [3J 2] , where the bipartite entangle- 
ment can be given a satisfactory quantitative definition 
in terms of entropy and its linearized versions, such as 
purity [5] [5] . The characterization of the global features 
of entanglement is more involved, unveiling in general 
different features of the many-body wave function [71 [8] , 
and it is becoming clear that the multipartite entangle- 
ment of a large system cannot be fully characterized in 
terms of a single (or a few) measure (s). 

Entanglement measures the nonclassical correlations 
between the different components of a quantum system 
and unearths different characteristics of its many-body 
wave-function. When the quantum system is large, it be- 
comes therefore interesting to scrutinize the features of 
the distribution of some bipartite entanglement measure, 
such as the purity or the Von Neumann entropy. Besides 
being of interest for applications, this is an interesting 
problem in statistical mechanics. In [5] we tackled this 
problem by studying a random matrix model that de- 
scribes the statistical properties of the eigenvalues of the 
reduced density matrix of a subsystem A of dimension 
N (the complementary subsystem B having the same di- 
mension as A). In the limit of large system dimension 
N, we introduced a partition function for the canoni- 
cal ensemble as a function of a fictitious temperature. 
The role of energy is played by purity: different temper- 



atures correspond to different entanglement. The most 
important result of our analysis was the proof that the 
different regions of entanglement, corresponding to dif- 
ferent ranges of the fictitious temperature, are separated 
by phase transitions. 

One puzzle was left open in that paper: in the region 
of negative temperatures our solution suddenly ceased to 
exist at a critical /3 g where the average purity of subsytem 
A was ttab = 9/4 N (a phenomenon quite common in 
random matrix theory, as this critical point corresponds 
to tesselations of random surfaces, or 2-D quantum grav- 
ity). As the partition function exists for all /3 € K, the 
region of factorizable states, where ttab = O (1), was not 
covered. 

We will show in this paper that the solution in [9] be- 
comes metastable for any j3 < (in the scaling of [5]) 
and a new stable solution appears which interpolates 
smoothly from ttab = 2/A to ttab = 1 as j3 goes from 
to — co. Moreover, we will also study the metastable 
solution that is born at f3 = and follow it through the 
region > (3 > — oo. 

This paper is organized as follows. In Sec. [TT| we in- 
troduce the notation and set the bases of the statisti- 
cal mechanical approach to the problem. In Sec. |III| wc 
study positive temperatures, where at very low temper- 
atures we find very entangled states. Negative tempera- 



tures are investigated in Sec. IV where it is shown that 
two branches exist, a stable one associated to a par- 
tial factorization of the many-body wave function, and 
a metastable one which contains the 2D-quantum grav- 
ity point. The finite size corrections are investigated in 
Sec. [V] We discuss the implications of our results for 
quantum information in Sec.|VI|and we conclude in Sec. 
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VII by summarizing our findings and discussing them in 
terms of future perspectives. 



II. A STATISTICAL APPROACH TO THE 
STUDY OF BIPARTITE ENTANGLEMENT 

We start off by describing a statistical approach to the 
study of bipartite entanglement for large quantum sys- 
tems. We will tackle this problem by introducing a par- 
tition function [9]. 

Consider a bipartite system, composed of two subsys- 
tems A and B. The total system lives in the tensor prod- 
uct Hilbert space H = Ha <8> with dim Ha = N < 
dim "Ha = M. We assume that the system is in a pure 
state £ H. The reduced density matrix of subsystem 
A reads 



p A = Tt B \^){^\ 



(1) 



and is a Hermitian, positive, unit-trace N x N matrix. 
A good measure of the entanglement between the two 
subsystems is given by the purity 



N 



(2) 



whose minimum is attained when all the eigenvalues Xj 
are equal to 1/N (completely mixed state and maximal 
entanglement between the two bipartitions), while its 
maximum (attained when one eigenvalue is 1 and all oth- 
ers are 0) detects a factorized state (no entanglement). 

In order to study the statistics of bipartite entangle- 
ment for pure quantum states, we consider typical vector 
states \ip) [TD1 [TT], sampled according to the unique, uni- 
tarily invariant measure. The significance of this mea- 
sure can be understood in the following way: consider 
a reference state vector \ipo) and a unitary transforma- 
tion l?/^) = U\ipo)- In the least set of assumptions on U, 
the measure is chosen in a unique way, being the only 
left- and right-invariant Haar (probability) measure of 
the unitary group U(N 2 ). The final state will hence 
be distributed according to the measure mentioned above 
(independently of \ipo))- Notice the analogy with the 
maximum entropy argument in classical statistical me- 
chanics. By tracing over subsystem J5, this measure in- 
duces the following measure over the space of Hermitian, 
positive matrices of unit trace PHI fTTj 



dfj,(pj 



Vp A (det p A 



,M-N 



5(1 - Tt Pa ) 



where Afc are the (positive) eigenvalues of pA (Schmidt 
coefficients), we dropped the volume of the U(N) group 
(which is irrelevant for our purposes) and rjN = M — N 
is the difference between the dimensions of the Hilbert 



spaces Ha and Hb- In order to study the statistical be- 
havior of a large bipartite quantum system we introduce 
a partition function from which all the thermodynamic 
quantities, for example the entropy or the free energy, 
can be computed: 



Z A b = J dp(p A ) exp (-(3N a TT A B) ■ 



(4) 



where a is a positive integer (either 2 or 3, as we shall 
see) and (3 a fictitious temperature "selecting" different 
regions of entanglement. The value of a needs to be 
chosen in order to yield the correct thermodynamic limit 



N a (n AB ) =0{N 2 ), 



(5) 



since N 2 is the number of degrees of freedom of the ma- 
trix pa- Around the maximally entangled states (for j3 > 
[9]) we have (ttab) = O (1/N) so a — 3, while around 
separable states (for j3 < 0) we have (ttab) =0(1) and 
hence a = 2. In the following we will assume rj = 0, since 
this does not change the qualitative picture (the exten- 
sion to r] 7^ being straightforward but computationally 
cumbersome). 

Since the purity depends only on the eigenvalues of pa 
the partition function reads 



J AB 



f ^An(A*-Aj) a f (i-iS) 



-PN" J2i x i 



(6) 



which by introducing a Lagrange multiplier for the delta 
function yields 



Z AB = N 2 



d N X 



-oo 2tt JA;>0 

A physical interpretation of the exponent in the inte- 
grand of the partition function can be given as follows 
[L2"] : the eigenvalues of pa can be interpreted as a gas of 
interacting point charges (Coulomb gas) at positions Ai's, 
on the positive half-line and with a quadratic potential. 
The solution of these integrals is known (as Selberg's in- 
tegral) for the case in which the integration limits are 
-oo < Xi < +oo QU- 

The constraint of the positivity of the eigenvalues 
makes the computation of this integral far more compli- 
cated. Although a exact solution for finite N is unlikely 
to be found [53] (but see 03] for the first few mo- 
ments), the problem arising from the constraint on the 
positivity of the eigenvalues can be overcome in the large 
TV limit, as we will look for the stationary point of the 
exponent with respect to both the A's and £. In particu- 
lar, the contour of integration for £ lies on the real axis, 
but we will soon see that the saddle point for £ lies on the 
imaginary £ axis. It is then understood that the contour 
needs to be deformed to pass by this point parallel to the 
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line of steepest descent. For the saddle point we need to 
find the minimum of the free energy: 



/3F = pN a A ?~ 2 E ln \\-^HN 2 d 1 - ]T *i ) 

i i<j \ if 



starting from f3 = +00 to /3 = —00 as 



By varying F we find the saddle point equations 



(8) 



2/3A Q A 2 + 2V-^— -iN 2 £ = 0, (9) 

— Aj — A 



E A > = !■ ( 10 ) 



In the following sections we will separately analyze the 
range of positive and negative temperatures, and unveil 
the presence of two phase transitions for the system, a 
second order one at a positive critical (3 and a first order 
one at a negative critical /3. 



A. The global picture 



Before proceeding to a formal analysis of the phase 
transitions, it is convenient to give a qualitative picture 
of the behavior of the Schmidt coefficients as the tem- 
perature is changed. As the inverse temperature is de- 
creased, the density matrix of subsystem A changes as 
follows. As /3 -> +00 all eigenvalues are = 1/N (max- 
imally mixed state). As f3 decreases, we encounter two 
phase transitions: one at a positive critical value /3+ and 
one at a negative critical value /3_, both critical values 
being to be determined. For /?>/?+ all eigenvalues re- 
main (1/N), their distribution being characterized, as 
we shall see, by the Wigner semicircle law. After the first 
phase transition, for f3 + > ft > ft_, the eigenvalues, all al- 
ways (1/N), follow the Wishart distribution, divergent 
at the origin. Finally, after the second phase transition, 
for ft < j8_, one eigenvalue becomes 0(1), "evaporating" 
from the "sea" of eigenvalues 0(1/N): this is a signa- 
ture of the emergence of factorization in the many-body 
wave function, the eigenvalue (1) being associated with 
a significant separability between subsystems A and B. 
For ft — > — 00 the many body wavefunction is fully fac- 
torized. Pictorially, the typical eigenvalues vector evolves 



1 1 1 

N' N''"' N 
*** 

,3-s-+oo 



o 



.,0 



N 



+oo>/3>/3 + 



,0,...,0 



0(1), 



/3+>0>l3- 
1 

N 



,0,...,0 



,3„>^>-oo 

(1,0,. ..,0), 

v v ' 



(11) 



where the zeroes in the second and third lines mean an 
accumulation of points around the origin, and we will 
find that [in the scaling of /3 given by a — 3 in Eq. jHJ], 
P+ = 2 and P_ = -2.455/JV. 



III. POSITIVE TEMPERATURES: TOWARDS 
MAXIMALLY ENTANGLED STATES 

In this section we will consider the range of positive 
temperatures: < j3 < +00; in particular we will study 
the occurrence of a second order phase transition at j3 = 2 
[5]. We will use a novel more general method, that will 
enable us to find all solutions and can be easily extended 
to negative temperatures. 

From the expression of the partition function one can 
easily infer that when ft — > +00 the typical states belong- 
ing to this distribution are maximally entangled states 
and correspond to the case = 1/N, Vi € {1, . . . N}. It 
then follows that for this range of temperatures the right 
scaling exponent in is a = 3. 

In order to estimate the thermodynamic quantities of 
the system we solve the saddle point equations (|9])-( 10 ) in 
the continuous limit, by introducing the natural scaling 



A, = ^A(t 4 ), 



< U = - < 1. (12) 



In the limit N — > 00, Eq. ^ becomes 



(3X + P dX 
Jo 



A - A' 4=°' 



(13) 



which is a singular Fredholm equation of the first kind, 
known as Tricomi equation |14j . The function 



dt S(X - X(t)) 



(14) 



is the density of eigenvalues we want to determine. A 
similar equation, restricted at /? = 0, was studied by 
Page QI]. 
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FIG. 1: Solution domain for different values of temperature: 
for each value of ft (indicated) the relative full line encloses 
the region of the parameter space such that the eigenvalue 
density is positive. The line m — S corresponds to the positive 
eigenvalues condition. 



According to the Tricomi theorem [TTJ [T3] the solution 
of the integral equation (13) lies in a compact interval 



[a, b], (0 < a < b). Let us set 
a + b b — a 



< S < to. 



(15) 



We map the interval [a,b] into the interval [—1,1] by 
introducing the following change of variables: 



A = to + xS, 



(x) = p(X)6. 



We get 



with 



g(x) =-- + 85m + (3S 2 x 

TT \ 2 

whose normalized solution (J (f>dx = 1) is 
7T J_ 1 V 1 - x 2 y - x 



Try! — x 2 



By using the constraint (10), that is 
-i 

A 4>{x) dx = 1, 
we can fix the Lagrange multiplier to obtain 



(16) 



(17) 



(18) 



, (19) 



(20) 



1 



fiS 2 2(1 - to) 



x - (3S 2 x 2 



The physical solutions must have a density 4>(x) that is 
nonnegative for all x E (— 1, 1). Let us look at the points 
where the density vanishes, <j)(x) = 0. From (21 ) one gets 



1 / 1 — to 



x± 



p5 2 



(22) 
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FIG. 2: Contour plots of the free energy in regions of the 
parameter space such that 4>{x) > 0, for (a) P = 0, (b) (3 = 1 
and (c) /3 = 3. Darker regions have lower free energy. The 
arrows point at the minima. 



where 



1 — TO 



■f3S 2 



(36^ 
2 



(23) 



For /3 > one gets that A > for every to and S, and 
(f>(x) > for x € [£_,£+]. The level curves z-t = ±1 are 



5 



given by m = T 1 (<5, /3), where 

r±(5,/3) = i± | ( 1 



2 



(24) 



They are symmetric with respect to the line to = 1 and 
intersect at 5 = and at S = \/2//3. Therefore, the 
condition (—1, 1) C [xi, X2] implies that the points (5, m) 
should be restricted to a (possibly cut) "eye-shaped" do- 
main given by 



max 



{J, IT (<$,/?)} <m<r+(M), 



(25) 



(recall the constraint to > <5 in ( 15 ) that expresses the 
positivity of eigenvalues). The right corner of the eye is 
at 



(6, to) 



(26) 



and belongs to the boundary as long as f3 > 2. For (3 < 2 
the eye is cut by the line to = 5. See Fig. [T] 

Let us remark that all points inside the region cor- 
respond to solutions of the saddle point equations. In 
other words we have a two parameter continuous family 
of solutions. We will look at the eigenvalue density that 
minimizes the free energy density of the system. From 
Eqs. ^ and ( 10 1 with a — 3 by applying the scaling ( 12 ) 
we get 



F 



N ^ 



= u 



1<J 



s + lnN + O 



N 2 f3 ^ 



InN 



bx\X(ti) - X{tj)\ 



N 



f + lnN + O 



hxN 



Here. 



/= Jim lf N - \nN 



(27) 



(28) 



is the free energy density in the thermodynamic limit, 
which reads 



Pf = pu- s, 



(29) 



in terms of the internal energy density u and the entropy 
density s, 



\ 2 4>(x)dx, 



dx dy <j>{x)<t>(y)\n{5\x - y\). (30) 
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FIG. 3: Free energy on the boundary of the region of the do- 
main where <j>(x) > 0, for different temperatures (indicated). 
Dashed line: free energy / on the lower boundary of the eye- 
shaped domain; full line: free energy on the upper boundary. 
The sought minima of the free energy can be inferred from 
the graph and coincide with the dots in Figs. [T] and [2] 



In order to compute the entropy density one integrates 
the Tricomi equation (17) and obtains 



(j)(x)dx / (j)(y)dy\n\y - x\ 



dxcf>(x) ln(x + 1) 



-it I dx(f>(x) J g{y)dy. 



(31) 



We get 



ft 2 fl/i 4 

u(S,m,0) = l-(l-TO) 2 + y -^-, (32) 
,* m 2(1 -to) 2 /3 2 <5 4 , S , oo , 
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FIG. 4: (Up) Density of the eigenvalues for /3 > j3+ = 2. 
(Down) Density of the eigenvalues for /3 = 2, 4 and 10. In the 
temperatures range /3 £ [2, oo] the solution is given by the 
semicircle law. 



and thus 



pf(S,m,P) 



(3-/3(1 

8 2 S 4 
16 



2(1 



^ 2 



2 



-In 



(34) 



The contour plots of the free energy are shown in Fig. [2] 
Note that /, as well as u and s, is symmetric with respect 
to the line m = 1. This Z 2 symmetry will play a major 
role in the following. The only stationary point (a saddle 
point) of the free energy density / is at the right corner 
of the eye (26), see Figs, [I] and [2j Thus, the absolute 
minimum is on the boundary. 
For j3 > j3+, where 



/3+ = 2, 



(35) 



the point ( 26 ) is also the absolute minimum, whereas for 
< (3 < f3 + the absolute minimum is at the right upper 
corner of the allowed region, 5 = T^(S,f3), namely at 



6 3 S 

m = 5, with 8 1 1 

H 4 2 



0. 



(36) 



See the dots in Figs. [T] and [2] The behavior of the free 
energy at the boundaries of the allowed domain is shown 
in Fig. [3] for different temperatures. 

We will study the behavior of our system starting from 
high values of 8, that is low values of internal energy 
u (purity). The analysis of lower values of 8, down to 



8 = and even below, will be done in the next section. 



For 8 > 8+ , by setting from ( 26 1 m — 1 and 

P 



2 



(37) 



and recalling (21 ), one gets the semicircle law (see Fig. [4]) 

2 

7T 



4>{x) = -y/T 



(38) 



whence, by (16) 



P(A) 




— \/A — a\/b — A, 

7T 



b=l+5=l 




(39) 



(40) 



This distribution is displayed in Fig. [4] Observe that as 
8 becomes larger the distribution becomes increasingly 
peaked around 1. This simply means that all eigenvalues 
tend to 1/N in the natural scaling (12): for tempera- 



tures T = 1/8 close to zero the quantum state becomes 
maximally entangled. 

By plugging ( 26 ) into ( 32 ) and ( 33 ) we get for 8 > 8+ 



u = 1 



5 1 
1 



= 1 



1 

2/3' 



Inf = -\ -lln(2p), 
2 4 2 v ; 



and thus 



Pf - 4 + 3 - in 5 = 8+ 3 + - H28). 

HJ S 2 4 2 H 4 2 V ' 



(41) 
(42) 

(43) 



At higher temperatures < 8 < P+ the solution ac- 



quires a different physiognomy. By plugging ( 36 ) into 

pi) 



4>{x) 



2 

it 5 



1 



(l + (2-6)x), 



see Fig. [5] yielding, by ( 16 1 
4 



p(X) = 



Tib 2 



b- A 
A 



6-2 



2(4 - b) 



(44) 



(45) 



with b = 2d. This is a Wishart distribution. See Fig. [5] 
The change from semicircle to Wishart is accompanied 
by a phase transition (the first of a series!) as we shall 
presently see. 



The half width S = b/2 is related to 8 by ([36 ) 

4 2 

53 ~ £2' 



P = 



(46) 



which runs monotonically from P = 2 when 5 = 1 to 
8 = when 6 = 2. Moreover, it reaches a minimum 
equal to 



P 9 



2 

27 



(47) 
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FIG. 5: (Up) Density of the eigenvalues for ,9 = 1. (Down) 
Density of eigenvalues for /3 = 0, = 2/3, and /3 = /3 g = 
—2/27 (dashed). In the range of temperatures /3 G (/3 9 ,/3+), 
with /3+ = 2, the solution is given by the Wishart distribution. 
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FIG. 6: Plot of Eq. (461 for positive (solid line) and negative 
(dashed line) temperatures. The minimum j3 g = —2/27 is 
attained at S — 3. 



at S = 3. Therefore, the above solution can be smoothly 
extended down to g , which is slightly negative, but not 
below. See Fig. [6] We will study the solution for negative 
temperatures in the next section. Note, incidentally, that 
the inverse function of ( 46 1 can be explicitly written 



A 



(48) 



With A = {^Wg + y/l~ H, H B > 

For < 0+ the internal energy (average purity) u is 



0/0 a ) 1/3 - 



obtained by plugging ( 46 1 into (|32 



(49) 



Therefore, at = (<5 = 2) one gets u = 2, at 0+ = 2 
(S = 1) one gets u — 5/4, and at g = —2/27 (5 = 3) one 
gets u = 9/4 (see the next section for the significance of 
these values) . From ( 46 1 and ( 33 1 one can also compute 



the entropy and the free energy for < 



0f = 



9 5 
-4 + S 
9 9 



3 

S 2 



In 



11 , 5 

— I In - , 

5 4 2' 



(50) 
(51) 



in terms of the function 6(0) € (1,3] introduced in Eq. 



Notice that 0f is the generating function for the con- 
nected correlations of ttab- The radius of convergence 
in the expansion around = 0, namely 2/27, defines 
the behavior of the late terms in the cumulants se- 
ries. Another interesting observation is that the func- 
tion r(x) = u((3 — —x/2) is the generating function of 
the number of rooted non-separable planar maps with n 
edges on the sphere (Sloane's A000139 also in [P7 | f!8 ] ). 
namely 

r(x) = 2 + x + 2x 2 + 6x 3 + 22.x 4 + 91a; 5 + 408x 6 

+1, 938a: 7 + 9, 614x 8 + 49, 335a; 9 + ... . (52) 

The counting of rooted planar maps on higher genus sur- 
faces is an unsolved problem in combinatorics and we 
conjecture it to be related to 1/N corrections of our for- 
mulas. 

We are now ready to unveil the presence of the first 



critical point at 0^ 



2. Let consider the density of 



eigenvalues ( 38 ) and ( 44 1 (or their counterpart ( 39 1 and 



(45 1). The phase transition at + is due to the restora- 
tion of a 7L>2 symmetry P ("parity") present in Eqs. ( [32| , 
( 33 ) and (34 1, namely the reflection of the distribution 



p(X) around the center of its support (m = S = 6/2 for 
< 0+ and m = 1 for > 0+). For < + there are 
two solutions linked by this symmetry, and we picked the 
one with the lowest /; at + this two solutions coincide 



with the semicircle ( 39 1 , which is invariant under P and 
becomes the valid and stable solution for higher 0. In 
order to explicitly show the presence of a second order 
phase transition in the system for — + we look at the 
expression of the entropy density s = 0(u — /), which 
counts the number of states with a given value of the 
purity. The expression for < + is given in Eq. (50 1, 
while for > 0+ it is given in Eq. ( 42 ) . 

At — + we get 5 = 1 and s = -1/4 - In 2. On the 
other hand the first derivative of s with respect to S is 
discontinuous at 6 = 1. However, also as a function 



of S, as given by (46 1 and (37), has a discontinuous first 



derivative at S — 1. By recalling that 



ds 
d0 
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d0" 
dS 



(53) 




FIG. 7: Entropy and its second derivative with respect to /?. 
The entropy is continuous in /3+ while its second derivative 
presents a finite discontinuity. 



one easily obtains that the discontinuities compensate 
and in the critical region, /3 — > we have 



In 2 



t4 + «-M^, (54) 



where is the step function. The entropy s is continuous 
at the phase transition, together with its first deriva- 
tive, although the second derivative is discontinuous, as 



shown in Fig. [7j Notice that the entropy is unbounded 
from below when /3 — ¥ +oo. The interpretation of this re- 
sult is quite straightforward: the minimum value of ttab 
is reached on a submanifold (isomorphic to SU(N)/Zn 
[H?] ) of dimension TV 2 — 1 , as opposed to the typical vec- 
tors which form a manifold of dimension 2N 2 — N — 1 
in the Hilbert space %. Since this manifold has zero vol- 
ume in the original Hilbert space, the entropy, being the 
logarithm of this volume, diverges. 

Now we want to express the entropy density s as a 



function of the internal energy density u. From (49) and 
(41| we get 



1 



6* 
4 



¥ 



<5 < 1, 
1<8<2, 



(55) 



that can be easily inverted 



3 - V9 - Au, § < u < 2. 



(56) 



function of S 



From (50 1 and (42) one gets the entropy density as a 



In ; 



< 5 < 1, 



(57) 



f + f-Jr+lnf, 1<5<2. 



Finally, by plugging (56) into (57), we obtain the entropy 
of the submanifold of fixed purity, s 



{\nV)/N 2 



function of its internal energy u = Nttab- 



s(u) 



|ln(u-l) 



hi 



(i-y/i 



1<«< |, 



(58) 



, 2 ; 4 



< u < 2. 



This function is plotted in Fig. [HJ 

Let us discuss the significance of these results. The 
present section was devoted to the study of positive tem- 
peratures T = > 0. In this range of temperatures, 
the eigenvalues of the reduced density matrix of our 7V 2 - 
dimensional system are always of O (1/N). As a conse- 
quence, the value of energy (purity) in Eq. ^ 



N 



i 

N 



X 2 p(X)d\ = O 



(59) 



is always small: there is therefore a lot of entanglement 



in our system. There are however, important differences 
as purity changes (it is important to keep in mind that 
in the statistical mechanical approach pursued here, the 
Lagrange multiplier /3 fixes the value of energy/purity). 
When 1/N < ttab < 5/4N the eigenvalues are dis- 
tributed according to the semicircle law (Fig. [4]), while 
for 5/47V < ttab < 2/N they follow the Wishart dis- 
tribution (Fig. [5]), the two regimes being separated by 
a second-order phase transition. The value ~kab — 2/N 
corresponds to infinite temperatures /3 = and there- 
fore to typical vectors in the Hilbert space (according to 
the Haar measure). One is therefore tempted to extend 
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FIG. 8: Entropy density s versus internal energy density u 
N{ttab}- See Eq. p8l. 



these results to negative temperatures [9] and one can 
indeed do so up to ttab — 9/47V, corresponding to the 
slightly negative temperature g — —2/27. However, as 
we have seen, a mathematical difficulty emerges, as this 
value represents the radius of convergence of an expan- 
sion around = and no smooth continuation of this 
solution seems possible beyond g . In the next section 
we will see that two branches exist for negative 0: one 
containing the point — g and in which purity is al- 
ways of O (1/N) and one in which purity is of O (1). The 
latter becomes stable for sufficiently large — /3's through 
a first order phase transition. 

Before continuing, we remind that larger values of pu- 
rity, towards the regime ttab = 0(1) yield separable 
(factorized) states. We are therefore going to look at 
the behavior of our quantum system towards separabil- 
ity (regime of small entanglement). 



IV. NEGATIVE TEMPERATURES 

A. Metastable branch (quantum gravity) 

By analytic continuation, the solution at positive of 
the previous subsection can be turned into a solution for 
negative 0, satisfying the constraints of positivity and 
normalization. In this section we will study this ana- 
lytic continuation and its phase transitions, but we an- 
ticipate that this is metastable for sufficiently large —/3's 
(namely for < — 2.455/iV) and that it will play a sec- 
ondary role in the thermodynamics of our model. How- 
ever, our interest in it is spurred by one of its critical 
points, at /3 = —2/27 = g which corresponds to the 
so-called 2-D quantum gravity free energy (see [H]), pro- 
vided an appropriate double-scaling limit (jointly — > g 
and N — > oo) is performed. 

In more details, the eigenvalue density (45) at = 
2/27, i.e. <5 = 3 [see between Eqs. (pfe| and (|48l, 



a 

and Fig. [6] reads 



P(A) 



2 / (6 - A) 3 
277rV A : 



(60) 



and from (49 1 u = 9/4 (see Fig. [7]and|9|. The derivative 
at the right edge of eigenvalue density in Fig. [5] vanishes. 
By expanding ( 46 1 for 6 — > 3 







27 
16 



2 

81 

;(S- 



2187 

that is, by setting x 



{5 -if 
-3) 5 



16 (*-3) 3 + i^-3) 



729 

e>((5-3) 6 



729 



(61) 



^2(0 - /3 9 )/9 -> 0, 



. , . , 4 2 35 3 80 4 1001 5 , n . 6N 

S = 3 ( 1 +x + -x 2 + —x A + —x 4 + -^x° + (x b 



and therefore for 



18 



0f 



3 , 

log 

4 S 2 

8lV2 



27 



216 



(62) 



81 
16 



(0-0 9 f 



{0-0 g f' 2 + o((0-0 g y 



(63) 



In fact, if one relaxes the unit trace condition, our parti- 
tion function Z has been studied in the context of random 
matrix theories [16, before. The objects generated in this 
way correspond to chequered polygonations of surfaces. 
Our calculations show that the constraint Tr pa = 1 is 
irrelevant for the critical exponents in this region. 

However, this is not a real critical point of our Coulomb 
gas. As this is an analytic continuation of the solution 
obtained for > 0, we are not assured that this is indeed 
a stable branch. In the next section we will show that 
a first order phase transition occurs at a lower value of 
0, namely at ~ — 2.455/A^ in this scaling (and there- 
fore the exponent a needs to be lowered from 3 to 2 for 
negative 0). The new stable phase will take over for 
all negative 0, where = —00 corresponds to separa- 
ble states. However the analytic continuation described 
here, although metastable, exists for all negative and 
we can study in more detail the behavior of the eigenvalue 
density (21) and of its free energy (34). The solution is 



straightforward but lengthy and is given in the following. 
It is of interest in itself because, as we shall see, it en- 
tails a restoration of the Z2 symmetry that was broken at 
the phase transition at = 2 described in the previous 
section. 

Recall that the density <fi(x) must be nonnegative for 
all x € (—1,1). This condition for < gives x ^ 
(x+,x_), with x± given by (22)- (23). The level curves 
x± = ±1 are given by m 
while the level curves A = 
with 



Tf (ST0), with Tf in (|24|), 



are given by m 



r ± 
1 2 



r*(6,0) = l±5\ -0 1 + 



0S 2 



(64) 



They are symmetric with respect to the line m = 1 and 
intersect at 5 — and at S = \J—2/f3. Moreover, they 
are tangent to at the points 



(6, m) 



2 2 
— , 1 ± - 

30' 3 



2 

30 



(65) 
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FIG. 9: u — N(ttab) as a function of (inverse) temperature. 
Notice that u — 2 at fi — (typical states). In the fi — > 
oo limit we find the minimum u = 1. The critical points 
described in the text are at fi g = —2/27, u = 9/4 (left point) 
and fi+ — 2,u — 5/4 (right point). However this phase is 
unstable (dashed line) towards another phase as soon as fi < 
(in this scaling of fi). 
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FIG. 10: Metastable branch. Domain of existence for the 
solution (m, 5) for negative temperatures. 



FIG. 11: Metastable branch. Contour plots of the free energy 
in regions of the parameter space such that <j>(x) > 0, for (a) 
P = Pg, (b) fi — —1 and (c) fi = —5. Darker regions have 
lower free energy. 



as shown in Fig. 10 Therefore, the condition (—1,1) n 
[aJi,a?2] — ^ implies that the points (6, m) should be re- 
stricted to a (possibly cut) eye-shaped domain given by 



max {5, h-(/3,$)} <m< h + ((3,S), 



where 



'rf(s,p), o<s<^f- 



_2_ 

3/3 ' 



2_ 
3d' 



(66) 



(67) 



The right corner of the eye is given by 



(6, m) = 



(68) 



and belongs to the boundary as long as f3 < —2. For 
P > —4 the eye is cut by the line m = 8. See Fig. (TO] 



The contour plots of the free energy (|34|) are shown 
in Fig. 
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The free energy density / has no stationary 
points for /3 < 0. The behavior of the free energy at the 
boundaries of the allowed domain is shown in Fig. |12| for 
different temperatures. For /? < —2 the right corner of 



11 



the eye is the absolute minimum, whereas for —2 < (3 < 
the absolute minimum is at the right upper corner of the 
allowed region, namely at 



with 5 = h + (fi,5), 



(69) 



that is 

and 
5-1 = 



^ + --1 = 0, for-A</?<0, (70) 



t[l + /3*p), for-2</3<-^. (71) 



Note that ( 70 ) coincides with ( 36 ) and thus is the prolon- 



gation of the curve ( 46 1 which runs monotonically from 
/? = when 5 = 2 to its minimum j3 g = —2/27 at 5 = 3. 



On the other hand, ( 71 1 is given by the curves 



P = -J2 ± ^- S2 + 4S ~ 2 



1 

"8* 5 3 



^yJ(2 + V2-6)(6-2 + V2), (72) 



which run from (3 — f3 g when 5 — 3 (with derivative zero) 
up to /? = — 3/2 + y/2 when 5 = 2 + y/2 (with derivative 
-oo) and then from /? = -3/2 + y/2 when 5 = 2 + y/2 
(with derivative +oo) up to (3 = —2 when 5 = 1. See 
Fig.pl 



Let us look at the eigenvalue density (21 1. When f3 q < 



13 < the solution is obtained by plugging (36) into (21 ) 



7T0 V 1 + X s 



(73) 



with 2 < 5 < 3, and is Wishart. At j3 g one gets 5 = 3 
and 



2_ {1-xf 
37V 1 + x ' 



(74) 



whose derivative at the right edge x = 1 vanis hes. On 
the other hand, when -3/2 + y/2 < (3 < f3 g by Q one 
gets 



tf>(x) 



1 



ir5\/l — x 2 



-(d+ V-5 2 + 4J - 2) + 2(1 - 5)x V-S 2 +'15-2^ 



with 3 < 5 < 2 + y/2, while for -2 < (3 < -3/2 + y/2 

1 



tt5^/1 — a; 2 



i (<5 - +4^-2) + 2(1 - 5)x +(5+ v / S 2 +45-2) 



(75) 



(76) 



with 1 < 5 < 2 + y/2. Note that this eigenvalue density diverge both at the left edge x = — 1 and at the right edge 
x = +1. 

I 



At /3 = — 2 one obtains 5 = 1 and 



get 



2x 2 



TTy/l 



(77) 



One gets the above density for all f3 < —2, where the 
Z 2 symmetry is restored. The interesting behavior of the 
eigenvalue density as j3 is varied is displayed in Fig. |14| 

The values of (to, 5) [that define the eigenvalue domain, 
see Eq. ( 15 )] and the thermodynamic functions u (inter- 
nal energy density) and s (entropy density) are shown 
in Figs. 15 16] respectively. Their explicit expressions 



are given for positive temperatures in Sec. Ill while for 
negative temperatures are given in the following. 

In the gravity branch, for (3 g < (3 < (2 < 5 < 3) we 



= 5, 
= -5 
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4 
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5 ~ S 2 
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2 



In 



5 2 - ln 2- 



(78) 



Beyond a second order phase transition at the critical 
temperature g we get that for —3/2 + y/2 < (3 < (3 g 
(3 < 5 < 2 + y/2) both u and s increase together with 
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FIG. 12: Metastable branch. Free energy on the boundary 
of the region of the domain where 4>(x) > 0, for different 
temperatures (indicated). Dashed line: free energy f3f on 
the lower boundary of the eye-shaped domain; full line: free 
energy on the upper boundary. The sought minima of the 
free energy can be inferred from the graph and coincide with 
the dots in Figs. [To| and [TT] 
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FIG. 13: Metastable branch. /3 vs 6. See Eq. ( 72 1 



the eigenvalue density halfwidth 5, 



m 


s, 


P = - 
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25- 
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+ 48- 
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\/-5 2 + AS - 2, 
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In 



6 
2' 



\/-5 2 + 45 -2 

(79) 



and then decrease for -2 < f) < -3/2 + a/2 (1 < S < 
2 + \/2), 



/8/ 



*. /3 = ~-^>/-<S 2 + 4S-2 ) 
25 - |j 2 + |V-<5 2 +45-2, 

5 25 17 

2 ~ 45 + 852 ■ 
5 



-In 



(80) 



Finally, beyond another second order phase transition 
for f3 < — 2 (0 < 5 < 1), when the Z2 symmetry is 
restored, we get 



= 1, P = - 



5 2 ' 



u = 1 



5 2 = 1 - 



3 

2^5' 



In - - - = --In (-2/3) - -, 
2 4 2 y H> 4' 



Pf = 



|-ln^ = -^+/3+iln(-2/?).(81) 



Finally, we record the interesting behavior of the mini- 
mum eigenvalue a = m — 8: see Fig.[l7j For —2 < (3 < 2, 
a coincides with the origin (left border of the solution 
domain). This variable can be taken as an order pa- 
rameter for both the second order phase transitions at 
P = — 2 and at @+ = 2. The Z2 symmetry is broken for 
—2 < P < 2. Notice, however, that the gravity critical 
point at P g = —2/27 remains undetected by a. 

Let us briefly comment on the fact that the analytic 
continuation of the solution at positive P described in 
the previous subsection has not led us towards separable 
states. The eigenvalues remain of O (1/N) (and so does 
purity) even though the temperature can be (very) neg- 
ative (as P crosses 0). In order to find separable states 
we will have to look at the stable branch in the following 
subsection. 
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FIG. 14: Metastable branch. Eigenvalue density for (3 = -1/27, f3 g = -2/27, (3 < f3 g and /3 < /3_ = -2. From left to right, 
notice how the distribution (initially Wishart, whose derivative at the rigth edge of the domain diverges) gets first a vanishing 
derivative at the right edge, then develops a singularity there and eventually restores the Z2 symmetry that was broken at the 
phase transition at /? = 2 described in the previous section. 
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FIG. 15: Average and width of the solution domain m and S 
[Eq. ( 15 1] as a function of /3. Solid line: stable branch. Dotted 
line: metastable branch. 



B. Stable branch of separable states 

In this section we will search the stable solution of 
the system at negative temperatures. As anticipated in 
Sec. |nj from the definition Q of the partition function 
one expects that, for any N, as (3 — > — oo the system 
approaches the region of the phase space associated to 
separable states: here the purity is O (1) and the right 
scaling in Eqs. is a — 2. In other words, by 

adopting the scaling iV 2 for the exponent of the parti- 
tion function, we will explore the region f3 = 0(1/N) 
of the scaling TV 3 introduced for positive temperatures. 
Notice that the critical point /? = —2/27 for the solution 
at negative temperatures now reads (3 = —{2/21)N and 
escapes to — oo in the thermodynamic limit. 
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FIG. 16: Internal energy density u and entropy density s 
versus /3. Solid line: stable branch. Dotted line: metastable 
branch. 
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We will show that the solution ( 45 1 , according to which 



FIG. 17: Minimum eigenvalue a — m — S versus (3. Solid line: 
stable branch. Dotted line: metastable branch. 
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all the eigenvalues are 0(1/N), becomes metastable in 
the region of negative temperatures, and the distribu- 
tion of the eigenvalues minimizing the free energy is such 
that one eigenvalue is O (1): this solution in the limit 
f3 — » — oo will correspond to the case of separable states. 
By following an approach similar to that adopted for pos- 
itive temperatures, we will first look for the set of eigen- 
values {Ai, . . . , Xn} satisfying the saddle point equations 
([9])-(10l with a = 2, getting as in Sec. Ill a continuous 
family of solutions. We will select among them the set 
maximizing (f3 < 0) the free energy ([8]), with a = 2: 



iV 



In - A 2 - 



— y 

l<i<j<N 



In \ Xj — Xj 



(82) 



As emphasized at the beginning of this section, since we 
are approaching the limit (3 — > — oo the states occupying 
the largest volume in phase space are separable; we then 
define An = M as the maximum eigenvalue and conjecture 
it to be of order of unity, whereas the other eigenvalues 
are <D(1/N): 



X N =p = (1) . 



E 



1 



p. 



(83) 



1<«JV-1 



From this it follows that we need to introduce the natural 
scaling only for the first N — 1 eigenvalues in order to 
solve the saddle point equations in the continuous limit 
and then estimate the thermodynamic quantities: 



A, = (1- M ) 



X(U) 
N-l' 



(84) 



< U = 



N - 1 



< 1, Vi= 1,...,N- 1. 



In particular we will separately solve the saddle point 
equations j9j-(10l, corresponding to the minimization of 
the exponent of the partition function with respect to 
the first N — 1 eigenvalues and the Lagrange multiplier 
£, given, in the limit N — ¥ oo, by 



P 



p{x')dx' x 



A- A 



with 



r -i|(l-M)=0, 
Xp(X)dX = 1, 



(85) 



and we will then consider the condition deriving from the 
saddle point equation associated to p 

2p/3 + i£ = 0. (86) 

The function p introduced in ( [85] ) is the density of the 
eigenvalues associated to Ai,...Aat_i in (84 1 and has 



the same form ( 14 1 introduced for p(X) in the regime 
of positive temperatures. By the same change of vari- 
ables introduced in Sec. |III[ Eqs. (fTo])-([l6|, the solution 



of the integral equations ( 85 1 can be expressed in terms 
of <j)(x) = p(X)5: 



i 



1 - 



2x(m- 1) 



(87) 



and the Lagrange multiplier is £ = — i4(m— 1)/(<5 2 (1 — //)). 
The region of the parameter space (m, S) such that the 
density of eigenvalues <j> is nonnegative reads 



5, 1- 



< m < 1 



which is the same expression of the domain found for the 
range of positive temperatur es (|25| ) when f3 = 0, namely 
r±(<5,0) = 1 ± 5/2 (see Fig. [l^~which is the analog of 
Fig. [I]). This is consistent with the change in temperature 
scaling from N 3 in the case of positive temperatures to 
N 2 in the case of negative temperatures: we are "zoom- 
ing" into the region near j3 — > CP of the range of tem- 
peratures analyzed in [9] and Sec. Ill Summarizing, as 



could be expected from what we have shown for positive 
temperatures, the solution of the saddle point equations 
is a two parameter continuous family of solutions. We 
now have to determine the density of eigenvalues that 
maximizes the free energy of the system. From Eqs. (83) 
and ( 82 ) we get 



f N =p 2 ~ 



N 2 /3 



2J \n\Xj - Xi 



o 



l<i<j<N 



(89) 



and by applying the scaling (84) 
o 1 



N 



(3 



ln(l - p) + f Ic d(8,m,P) 



(90) 



where 



/= lim /jv--lniV 



ln(l-/i) + /rod, 



(91) 

is the free energy density in the thermodynamic limit, 
and 



f rcd (5,m,(3) = ~J dx(f>(x) 

2(m-l) 2 1 (5 



dy<f>(y) ln(S\x - y\) 



(92) 



is the reduced free energy density of the sea of eigenval- 
ues. 

It is easy to see that Pf Te d{m, <5), has no stationary 
points, but only a global minimum /3/ rc d = 1/2 at 
(S, m) = (2,2), see arrow in Fig. 18 this point yields 
the Wishart distribution found at f3 — for the case of 
positive temperature (see also [S]): 



1 1-x 



1 / 4 - A 
2tt" 



(93) 



7T V 1 + x 1 rv/ 2n V A 
where one should remember that the A's are also scaled 



by 1 - p, see Eq. (84). 
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FIG. 18: Contour plot of the reduced free energy (3f le d(8,m) 
of the sea for negative temperatures. The arrow points at the 
minimum. 
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FIG. 19: Reduced free energy /3/ ro d on the boundary of the 
triangular domain in Fig. [18] for the case of negative tem- 
peratures. Solid line: upper boundary; dashed line: lower 
boundary. 



We stress that this result is valid for all /3 < 0. In order 
to check this solution one has to compute the free energy 
on the boundary of this domain, see Fig. [l9| (which is the 
analog of Fig. [3]). One gets that the free energy density 
is given by 



f(t*,P)=P? 



P 



In(l-M) 



1 

w 



(94) 



A new stationary solution, in which the largest isolated 
eigenvalue /i becomes O (1), can be found by minimizing 
the free energy density and yields 



KP) 



i 



2 ' \t + 1' 



(95) 



being defined only for f3 < —2; this expression can also be 
obtained directly by the saddle point equation ( 86 ) cor- 



responding to the isolated eigenvalue fi. This eigenvalue, 
O (1), evaporates from the sea of eigenvalues O (1/N), as 



FIG. 20: Evaporation of the eigenvalue fj, = O (1) from the 
sea of eigenvalues O (1/N). 




FIG. 21: Reduced free energy as a function of for different 
values of /?(< 0). Notice the birth of a stationary point for 
/? = — 2 that becomes the global minimum for /3 < /?_. 



pictorially represented in Fig. [20] The isolated eigenvalue 
moves at a speed -dn/d/3 = l/(2 v / /3 4 + 2/3 3 ), which di- 
verges at f3 = —2: another symptom of criticality. How- 
ever, this new solution, when it appears at /3 = —2, is not 
the global minimum of /?/: as we shall see it eventually 
becomes stable at a lower value of /3. We get for /3 < — 2 
(i.e. < fi < 1) 



ln(l - M ) 
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= In 
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(3f = 



1 - 2/i 

P M — s = 777^ v — m (1 — /i) 
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(96) 



,(97) 



(98) 



We are now ready to unveil the presence of a first or- 
der phase transition in the system. In Fig. [21] we plot 
the free energy density as a function of /i for different 
values of j3. For j3 > — 2 there is a global minimum of fif 
at i.i = 0; /i is still in the sea of the eigenvalues O (1/N) 
and the stable solution is given by the Wishart distri- 
bution (73 1 with the potentials (78) (remember that, in 
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FIG. 22: Free energy and maximum eigenvalue at negative 
temperatures. The two solutions are exchanged at /3_ ~ 
2.45541, where there is a first order phase transition. Full 
line: solution of mimimal free energy; dashed line: solution of 
higher free energy. 



the zoomed scale considered here, j3 g corresponds to the 
very large inverse temperature Nj3 g ). At f3 — — 2 there 
appears a stationary point for the free energy density 
corresponding to = 0(1) [see (95)]; notice however 
that Pf at this point remains larger than its value at the 
global minimum, until P reaches /3_. Finally, for 
the global minimum of Pf moves to the right, to the so- 
lution containing \i — O (1). Summarizing, for j3 > /5L 
the solution of saddle point equations maximizing the 
free energy of the system is such that all eigenvalues are 
0(1/N), at P = —2 there appears a metastable solu- 
tion for the system with one eigenvalue 0(1), and for 
P < P- this becomes the stable solution, that maximizes 
the free energy, whereas the distribution of the eigen- 
values found in Sec. IIHI becomes now metastable. The 
maximum eigenvalue is then a discontinuous function of 
the temperature at /3 = /3_ and in the limit /3 — > — oo, 
/i approaches 1: the state becomes separable. This criti- 
cal temperature /3_ is the solution of the transcendental 
equation /(/?_, 0) = /(/3_,/x_), that is 



2(1 -M- 



which yields 
//_ ~ 0.71533, 



= -ln(l- M -), 



2M1-M-) 



(99) 



-2.45541. 



(100) 



0. 



hand, the solution [i = 0, corresponding to 
u = /i 2 = 

s = P(u 

M - \ 

has a lower free energy for / 



/) 



1 

2' 



(101) 
(102) 

(103) 

< P < 0, and a higher one 
for P <P_. See Fig. [22 

At P- there is a first order phase transition. At this 
fixed temperature the internal energy of the system goes 
from u r = up to ui = (jr. ~ 0.5117, while the entropy 
goes from s r = —1/2 down to sj = —1/2 + ln(l — /i_) ~ 
— 1.75643. One gets As/Au = Therefore, the en- 
tropy density as a function of the internal energy density 
reads 



s(u) 



P-U - |, 



ln(l - y/jj) 



< u < n 2 _ , 



(J£_ < U < 1. 



(104) 



It is continuous together with its first derivative at u = 
while its second derivative is discontinuous. Notice 
that Au = As IP- is the specific latent heat of the evapo- 
ration of the largest eigenvalue from the sea of the eigen- 
values, from (l/N) up to 

A few words of interpretation are necessary. As we 
have seen, it has been necessary to follow the stable 
branch of the solution in order to obtain separable states 
at negative temperatures. The analytic continuation 
of the stable solution for positive temperatures would 
yield an unstable branch in which all eigenvalues remain 
0(1/N). By contrast, the new stable solution consists 
in a sea of N — 1 eigenvalues (l/N) plus one isolated 
eigenvalue (1). 

Let us discuss this result in terms of purity, like at the 



end of Sec. Ill (we stress again that P is a Lagrange multi- 



plier that fixes the value of the purity of the reduced den- 
sity matrix of our iV 2 -dimensional system). Assume that 
we pick a given isopurity manifold in the original Hilbert 
space, defined by a given finite value ttab of purity. If 
we randomly select a vector belonging to this isopurity 
manifold, its reduced density matrix (for the fixed bi- 
partition) will have one finite eigenvalue \i ~ ^Jt^ab and 
many small eigenvalues 0(1/N) (yielding a correction 
(l/N) to purity). In this sense, the quantum state is 
largely separable. The probability of finding in the afore- 
mentioned manifold a vector whose reduced density ma- 
trix has, say, two (or more) finite eigenvalues /ii and /i2 
(such that + [i\ ~ ttab, modulo corrections O (l/N)) 
is vanishingly small. By contrast, remember (from the 
results of Sec. Ill ) that if the isopurity manifold is char- 



acterized by a very small value 0(l/N) of purity, the 
eigenvalues of a randomly chosen vector on the mani- 
fold are all (l/N) (being distributed according to the 
semicircle or Wishart, depending on the precise value of 
purity, as seen in Sec. III). This is the significance of the 



Therefore, the branch (97)- (98 1 is stable for P < P_ while statistical mechanical approach adopted in this article. 



it becomes metastable for P_ < P < —2. On the other We will come back to this point in Sec. VI 
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FIG. 23: Finite N version of Fig. |22| Free energy and maximum eigenvalue in the saddle point approximation as function of 
at N = 30. The local minimum is in blue, the global one in red. The two minima swap stability at = —1.935. Notice 
the birth of the new local minimum at — —1.8 (for N = oo this takes place at = —2) and the exchange of stability at 
= -1.93 (for N = oo, = -2.45). 



V. FINITE SIZE SYSTEMS 

The results of the previous section refer to N —> oo. 
In order to understand how finite-TV corrections affect 
our conclusions we have numerically minimized the free 
energy for various temperatures. The two phases of the 
system discussed in the previous section correspond to 
the two solutions obtained by minimizing the free energy 
PJn p2J| on the N dimensional simplex of the normal- 
ized eigenvalues. Indeed, we have numerically proved 
that PJn{P) presents two local minima at negative tem- 
peratures: for > j3 > B^ the minimum giving the 
lower value of PJn{P) corresponds to the distribution of 
eigenvalues ( 45 ) , found in the last section; the other min- 



imum is reached when the highest eigenvalue is 0(1). 
The point B — P_ is a crossing point for these two 
solutions, and for 8 < /3_ these two solutions are in- 



verted, see Fig. [22] (and 21 1. Summarizing, there exists 



a negative temperature at which the system undergoes 



a first-order phase transition, from typical to separable 
states. 



The first thing to notice is that qualitatively the phase 
transition remains of first order even for finite N. The 
second is that the finite N corrections are quite relevant 
for the location of the phase transition and the value of 
the maximum eigenvalue as a function of p. For example, 
for N — 30, the negative critical temperature p^ ^ — 
— 1.935 instead of —2.455. This is evinced from Fig. 23 
which is the finite size version of Fig. [22] This can be 
understood, as the corrections to /(/i) around /i = are 
quite large. In the limit \i = 1/N there is a hard wall for 
the maximum eigenvalue /i, as the condition J^. Aj = 1 
cannot be satisfied if /i < 1/N. It is therefore likely 
that all sorts of large corrections occur as [i tends to 
1/N, probably yielding an effective size to the corrections 
which is a lower power of 1/N (or even possibly 1/ In AT). 
The limits [i — > and TV — > oo do not commute. 



To further explore this effect we minimized /3/at with 
respect to Ai, . . . , \n-i, for fixed values of the largest 
eigenvalue Ajv = A* an d for different temperatures. The 
results for N = 30 are shown in Fig. [24] One can 
see that between f3 = —1.8 and j3 — —2.5 there is a 
competition between two well defined local minima, cor- 
responding to the two solutions discussed above. At 
P — /3^ 30 ' = —1.935 their free energies are equal. For 
highe r 8 the global minimum corresponds to the solu- 
tion (45), whereas on the other side of the solution 
with /i = O (1) minimizes /3/jv- Similar corrections are 
observed for the value of 8 = at which the second 



minimum is born. See Fig. |25| 

We have seen that for 8 > j3- the stable solution has no 
detached eigenvalue. By taking into account the scaling 
8 — > B/N we get that the solution is given by the very 
first part of the gravity branch (78). In particular, the 
maximum eigenvalue is given by b/N — (m + S)/N = 
25 /N. On the other hand for B < 8_ the maximum 



eigenvalue is given by (95). Therefore, we get 



0<P-, 
P- < 8 < 0, 



(105) 
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'0.0 0.2 0.4 0.6 0.8 1.0 
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FIG. 24: Finite TV version of Fig. |21| /3/jv — In N as a, function of the maximum eigenvalue fj,, obtained by numerical minimization 
over the remaining TV — 1 eigenvalues for various f3. Observe the formation of a new minimum and the exchange of stability, 
although the critical values of /3 at which these phenomena occur differ from the theoretical ones, due to large finite TV 
corrections. However, it is clear that at small fi, 1/TV corrections tend to increase the value of /3/jv, making the critical value 
P- move towards 0, as observed in the numerics. 
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■1.75 



■U 



1.85 



25 



30 



35 



N 



with 8(f3) given by (46)-(48). The numerical results for 
N = 40 are compared with the expressions in Eq. ( 105 ) 
in Fig. [26] The agreement is excellent. 



FIG. 25: Value of /3 at which the second minimum is born 
as a function of TV. The solid curve is a best fit returning 



(N) 



-1.997 - 6.04/TV. The asymptotic value should be 2 



and is in good agreement with the constant of the fit. 



The corresponding free energy follows from ( 98 1 and 



( 78 1 with the appropriate scaling 



+ f+V 1+ i- ln (^v 1+ t)' / ^ /3 -< 

Pf = { (106) 

4 6{P/N) + S(/3/N) 2 111 2 ' p_ ^ p ^ U. 



Notice that in order to have a finite size scaling of the 
critical temperature one should take into account 

0(1/N) corrections to the expression of f3f and then 



evaluate the intersection between the two branches, but 
this analysis goes beyond our scope. 
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Ntt AB 

FIG. 27: Entropy density s versus internal energy density u. Notice that the unit on the abscissae is 1/N in the left panel. 



iV = 40 
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FIG. 26: Maximal eigenvalue in the saddle point approxima- 
tion as function of /3. The points are the result of a numerical 
evaluation for N = 40, while the full line is the expression in 



Eq. (1051 



VI. OVERVIEW 

Let us summarize the main results obtained in this ar- 
ticle in more intuitive terms, by focusing on those quan- 
tities that are more directly related to physical intuition. 



In the statistical mechanical approach adopted in this ar- 
ticle, the temperature plays the usual role of a Lagrange 
multiplier, whose only task is to fix the value of energy 
(purity in our case). A given value of /3 determines a set 
of vectors in the projective Hilbert space whose reduced 
density matrices have a given purity (isopurity manifold 
of quantum states). The distribution of the eigenvalues 
of (the reduced density matrices associated to) these vec- 
tors is that investigated in this article and yields infor- 
mation on the separability (entanglement) of these quan- 
tum states. The distribution of eigenvalues is the most 
probable one [T^] (in the same way as the Maxwell dis- 
tribution of molecular velocities is the most probable one 
at a given temperature). Let us therefore abandon tem- 
perature in the following and fully adopt purity as our 
physical quantity. 

Entropy counts the number of states with a given value 
of purity and is in this sense proportional to the loga- 
rithm of the volume in the projective Hilbert space. The 
explicit expressions of the entropy density s, which is the 
logarithm of the volume of the isopurity manifold, as a 
function of the purity ttab of the state vectors in that 
volume, can be read directly from Eqs. ( 58 1 and ( 104 1 by 
taking into account the correct scaling: 



s(ttab) 



| HNn AB - 1) - 1 
In (| - yf\ Nn AE 



fi-KAB — 2> 

.In (l - ^ir AB ) - §, 



2(IVPH i(i-Vi-^A B ) 



4N ' 



2 : 4JV 



< T^AB < 



N ' 



2 < TTAB < M- 



,2 



fit < 1T AB < 1, 



(107) 



with /£ ~ 0.512 and /3_ ~ -2.455 given by (99 )-( 100 1. The plot of s vs ir AB in the two regions ir AB = O {1/N) 
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and ttab = O (1 ) is shown in Fig. 27 By exponentiating 
the expression (107) we get the volume V = exp(iV 2 s) 



(i.e. the probability) ol the isopurity manifolds 



e-^r(Nn AB -l) N2 /\ 



1 ^ <f 5 



N 



AN ' 



V{-k A b) OC < 



(t-^f-A^) 

exp [iV 2 (/3— 7TAB -|)] 

ye 2 (1 - ^J-kab) , 



exp 



iV 2 - 



4 2(1-^1-^1) 4(|-V!-JV7r AB )' 



< < 



4j¥ 



(108) 



AT 



ln(VQ " 

-2 
-3 




ln(^) 
In(TTAB) 

FIG. 28: Volume V = exp(iV 2 s) of the isopurity manifolds 
versus their purity ttab for N = 50. 




1 1.1 5/4 1.4 NlT AB 

FIG. 29: Minimum eigenvalues as a function of ttab- At 
ttab = 5/4iV the gap vanishes. 

This is plotted in Fig. [28] for N = 50. The presence 
of discontinuities in some derivatives of entropy detects 
the two phase transitions. At ttab — 5 /AN there is a 
second order phase transition signaled by a discontinuity 
in the third derivative. Indeed, in general, if s, u and T 
are entropy, energy and temperature, respectively, and 




FIG. 30: Maximum eigenvalue as a function of ttab- 



C = du/dT is the specific heat, one gets ds/du = f3 = 
1/T and d 2 s/du 2 = -1/(T 2 C) = -(l/T^ds/dT)- 1 . 
Discontinuities of the nth derivative of ds/dT translate 
therefore in discontinuities of the (n + l)-th derivative 
of ds/du. The first order phase transition, which takes 
place between ttab = and ttab — M- — 0.512 is 
signaled by discontinuities in the second derivative of the 
entropy at those points. Observe that entropy is un- 
bounded from below: at both endpoints of the range of 
purity, tt ab = 1/N (maximally entangled states) and 
ttab = 1 (separable states), when the isopurity mani- 
fold shrinks to a vanishing volume in the original Hilbert 
space, the entropy, being the logarithm of this volume, 
diverges, and the number of vector states goes to zero 
(compared to the number of typical vector states). See 
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FIG. 31: Summary (stable branch): profiles of the eigenvalue density for 1 < ttab < 



< ttab < 



and 



< TTAB < 1. 



Fig. [28} 

The presence of the phase transitions can be easily read 
out from the behavior of the distribution of the Schmidt 
coefficients, i.e. the eigenvalues of the reduced density 
matrix pa of one subsystem. From Eq. (56) we get the 



expression of the minimum eigenvalue A m ; n = a = m — 5 
as a function of ttab 



Amin 



' £ (1 - 2y/NTTAB ~ 1) , jj<TTAB<ik, 



0, 



477 < 7TAB < I- 



(109) 



which is shown in Fig. |29| The second order phase transi- 
tion at ttab = 5/4iV, associated to a Z 2 symmetry break- 
ing, is detected by a vanishing gap. On the other hand, 
the maximum eigenvalue A max coincides with the upper 



edge of the sea of eigenvalues b = m + 5, as given by (56 1, 
until it evaporates according to Eq. (96). Thus, 



f i (1 + 2y/NTT AB ~ 1) , 4 < TTAB < 



Amax = S ^ (3 — 2W J — N-KAB ) , 47v 



5 

47V ' 



5 < *AB < 2 



AT 



< n AB < 1) 



(110) 



as shown in Fig. [30] 

In the different phases the distribution of the eigen- 
values of pa have very different profiles. See Fig. [3l"| 
While for 1 < ttab < 477 the eigenvalues (all O (1/N)) 
follow Wigner's semicircle law, they become distributed 
according to Wishart for larger purities, ^ < ttab < jj, 
across the second order phase transition. This is a first 
signature of separability: some eigenvalues vanish and 
the Schmidt rank decreases. For even larger values of 
purity, < ttab < 1) across the first order phase tran- 
sition, one eigenvalue evaporates, leaving the sea of the 
other eigenvalues O (1/N) and becoming 0(1). This is 
the signature of factorization, fully attained when the 
eigenvalue becomes 1 at ttab = 1- 

We tried to give an overview of the phenomenology of 
these phase transitions in Fig. [32j where we also showed 
the presence of the metastable branches discussed in 
Sec. |IV| The global picture is both rich and involved and 
it would not be surprising if additional features would be 
unveiled by future investigation. 



VII. CONCLUSIONS 

We have obtained a complete characterization of the 
statistical features of the bipartite entanglement of a 
large quantum system in a pure state. The global picture 
is interesting as several locally stable solutions exchange 
stabilities. On the stable branch (solutions of minimal 
free energy) there is a second order phase transition, as- 
sociated to a Z 2 symmetry breaking, and related to the 
vanishing of some Schmidt coefficients (eigenvalues of the 
reduced density matrix of one subsystem), followed by a 
first order phase transition, associated to the evaporation 
of the largest eigenvalue from the sea of the others. 

In the different phases the distribution of the Schmidt 
coefficients have very different profiles. While for large 
f3 (small purity) the eigenvalues (all 0(1/N)) follow 
Wigner's semicircle law, they become distributed accord- 
ing to Wishart for smaller (3 and larger purity, across the 
first transition. For even smaller (and eventually neg- 
ative) values of /3, when purity becomes finite, across 
the second phase transition, one eigenvalue evaporates, 
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FIG. 32: Overview of the "evolution" of the eigenvalue densities as a function of purity (1/N < ttab < 1)- The straight 
segment represents the stable branch: the distribution starts as a delta function, evolves into a semicircle, undergoes a second- 
order phase transition at ttab = 5/4jV (Z2 symmetry breaking), becomes Wishart and undergoes a first-order phase transition 
between 2/N < ttab < 0.5117, during which one eigenvalue evaporates from the sea of the others O (1/N) and becomes O (1). 
A metastable branch is born at ttab = 2/N: it starts as Wishart, undergoes a second-order phase transition at ttab = 9/4N 
(2-D gravity), develops a singularity at its right edge through a second order phase transition at ttab ~ 9/4N, then its support 
starts decreasing, undergoes a second-order phase transition at ttab = 7/4N (Z% symmetry restoration) and eventually becomes 
sharply peaked (with two singularities). The diamonds indicate the three second order phase transitions. 



leaving the sea of the other eigenvalues {1/N) and be- 
coming 0(1). This is the signature of separability, this 
eigenvalue being associated with the emergence of factor- 
ization in the wave function (given the bipartition). This 
intepretation is suggestive and hints at a profound mod- 
ification of the distribution of the eigenvalues as f3, and 
therefore purity, are changed. Remember that /3, viewed 
as a Lagrange multiplier in this statistical mechanical ap- 
proach, localizes the measure on set of states with a given 
entanglement (isopurity manifolds [19]). 

It would be of great interest to understand whether the 
phase transitions survive even in the multipartite entan- 
glement scenario, if one views the distribution of purity 
(over all balanced bipartitions) as a characterization of 
the global entanglement of the many-body wave function 
of the quantum system [5] . This description of multipar- 
tite entanglement displays the symptoms of frustration 
[20] . catapulting the problem into one of the most fasci- 



nating arenas of modern statistical mechanics [2TJ . 

While we were completing this work a paper appeared 
in which an aspect of this problem is discussed, although 
with a different emphasis [25]. In order to connect our 
results to those in [22] , notice that the probability distri- 
bution of the purity P(u) is proportional to the volume 
of the isopurity manifold and therefore 



lnP(u) ~ N 2 s(u), 

where s is the entropy and N(ttab) 
energy. See Eq. (|108[). 



(Ill) 



u is the internal 
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